Enhancement of condensate depletion due to spin-orbit coupling 
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We show that spin-orbit couphng(SOC) significantly enhances the depletion of a homogeneous 
Bose-Einstein condensate in three dimension. With decreasing anisotropy of SOC, both quantum 
and thermal depletion increase. Particularly, different types of SOC give rise to qualitatively different 
dependences of condensate depletion on microscopic variables including scattering length, strength 
of SOC and temperature, a novel feature that can be directly observed once these types of SOC are 
realized in experiments. Moreover, we point out that thermal depletion in three dimension becomes 
logarithmically divergent at any given finite temperature when both SOC and interaction approach 
the isotropic limit. 

PACS numbers: 



The recent realization of synthetic spin-orbit cou- 
pling(SOC) for ultra cold atoms is an exciting develop- 
ment in the field of quantum gases[l-4 . As many ex- 
perimental parameters, including density, interaction and 
the configuration of SOC itself [5 , are highly controllable, 
quantum gases provides physicists an ideal platform to 
investigate the interplay between SOC and interaction in 
many-body systems, a challenging problem that remains 
unsolved in condensate matter physics so far 

In the literature, extensive studies have predicted a 
number of novel properties of Bose-Einstein condensate 
in the presence of SOCfT'-TF . However, condensate de- 
pletion, an intrinsic property of bosons, has not been sys- 
tematically discussed except for a few numerical studies 
for Rashba SOC[l4l|T5]. In the absence of SOC, one im- 
portant effect of condensate depletion is that it destroys 
a condensate at low dimensions. In three dimension, 
the depletion for weakly interacting bosons is small and 
a condensate naturally exists at low temperatures [17 . 
However, in the presence of SOC, due to the lack of a 
systematical study of analytic expressions for the conden- 
sate depletion, it is unclear how the depletion depends on 
microscopic parameters, and in particular, on the config- 
uration of SOC. In this Rapid Communication, we point 
out that condensate depletion becomes significant even 
in three dimension in the presence of SOC and may com- 
pletely destroy the condensate at low temperatures. We 
also obtain the explicit forms for the dependence of con- 
densate depletion on the anisotropy and strength of SOC, 
the scattering length, and the temperature. 

The single-particle Hamiltonian with a general SOC 
between two spin states | t) and | 1) can be written as 

^o = ^{5I(-^a-2i«aa„a„) + 4}. (1) 

a 

where > i^y^i^z has been assumed without loss of 
generality. Throughout the paper we set h = 1. tZa 
is the strength of SOC, and crQ,(a = x^y^z) are Pauli 



matrices. The eigen-state has two branches |k^) = 
^/k|kt) +^ke'^^|k;), |k-) = -^ke-^^^|kt) +^/k|k;), cor- 
responding to eigen-energy 

e^ = ^{kl+kl + kl±2A^,^ + nl), (2) 

with A^,k = ^Klkl + Klkl + Klkl, Uk = ^^5(1 + ^), 

vvi = ^"^"^ "^k = arg{K^kx + iKyky). 

For isotropic SOC where 

— — — ^1 Eq. (j2j) 

becomes = 2^(|k| ± /^)^. The lower branch (e^) 
gives rise to the single particle density of states (DOS) 
Po(e) ^ ^ ~^ 0[18l [19], the same as that for 

one dimensional systems in the absence of SOC. As con- 
densate depletion of non-interacting bosons solely relies 
on single-particle DOS, we conclude a non-interacting 
condensate in three dimension is completely destroyed by 
the isotropic SOC even at zero temperature. The same 
conclusion applies to Rashba coupling in two dimension, 
where k,x = f^y 

While the above conclusion on condensate depletion 
of non-interacting bosons is unambiguous, it is essential 
to ask what are the condensate depletion in the pres- 
ence of both SOC and interaction, and how the depletion 
depends on the anisotropy of SOC, which always exists 
in practice. We consider the interaction between spin- 
orbit coupled bosons in the mean- field region [20], Hint — 
J d^r{2U^ih^{r)hi{r) -\-U^^ ^li^))^ where n^(r) is the 
density operator for spin cr(=t, i) and Uaa' is the interac- 
tion strength. For weak interactions, chemical potential 
is much smaller than the energy difference between the 
bottoms of two branches, i.e., ji <^ = Ai:^/2m, the 
negligible occupation in the higher branch justifies that 
Hint can be expanded in the lower branch (aj^,ak), as 

TJ ^ ^k3,k4 ft x 

TLint - ^ 2^ /ki,k2%i^k2^k3<^k4^)ki+k2-k3-k4, 

ki,k2,k3,k4 

(3) 
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FIG. 1: (Color Online) DOS and condensate depletion for three particular configurations of SOC. Red spheres represent the 
condensates, and orange clouds represent the depletion, which is enhanced along the y direction in (b) and the y and z directions 
in (c). Purple circle in (b) and sphere in (c) represent the kinetic energy minimums. Short blue arrows represent the sound 
modes and long green arrows represent the soft modes. The effective dimension (Deff) is obtained by comparing DOS with 
that for a Deff dimensional system without SOC. 



where 



+ ^^ki^k2^k3^k4 



o*(0ki +0k2 -0k3 -0k4) 



(4) 



U 



U has been assumed 



and Q is the volume. U^^ 
to simplify notations. 

Bogoliubov spectrum It has been shown that the 
mean field ground state in two dimension can be a plane- 
wave condensate or a stripe phase [7 . In three dimension, 
the situation is similar but more sophisticated [2T. In 
this Rapid Communication, we concrete on the plane- 
wave condensate at ko = k.x^x^ the mean- field ground 
state when U > U^^ and Kx ^ 

We perform Bogoliubov theory to study the quasi- 
particle spectrum, which can be written as £^±,q = 
± (e+,q - e_,q)/2, where 



)V4- 



(5) 



and e±,q = eko±q + f (/7 + 2/7^;iiko±q^ko±q cos 0ko±q), 
Vq = f [Wo+q^ko-qe'^^^O+c, + 0ko-.) + [/^Xko+q^ko-q + 



t^n(^ko+q^ko-qe'^^0+^ + ^ko+q^ko-q^ 



0]. For 



discussions, it is useful to enumerate £^q for three par- 



ticular configurations of SOC, (a) i^x 7^ 0, f^y 



0, 



(b) i^x i^y, i^z 0, (c) i^x = i^y = i^z, U U^i. (a) is 
the case when SOC exists along only one direction, (b) 
is the Rashba coupling, and (c) is referred to as isotropic 
SOC and interaction. For these three cases, £^q in low 
momentum limit have distinct expressions, 

(a) : 

(b) : 













(6) 



where q = |q|, the chemical potential ji = nil ^ and U = 
{U + U^i)/2. Since e+,q — e_,q ~ o{q^)^ E^^ represents 
the quasi-particle spectra in the long wavelength limit. 
For case (a), £^q is linear along all three directions and 
characterized by an isotropic sound velocity, i.e., E^^ = 
cq^ where c = \/fi/m the same as the sound velocity 
for an ordinary condensate without SOC. For Rashba 
coupling in case (b), the dispersion becomes quadratic 
along the y direction, i.e., a soft mode emerges(see also 
[ini HH US] ) • In case (c), the second soft mode emerges 
along the z direction, and this requires that both SOC 
and interaction are isotropic These soft modes have 
profound effects on condensate depletion. 

Depletion and DOS of Quasi-particles In the 
framework of Bogoliubov theory, condensate depletion is 
given by n^^ = ng^ 



where 



-Yi 



2n 



2E„ 



1), 



(7) 



-y 



^+,q ■ 



-,q 



2^n 



1 



^0 



and 



are quantum and thermal depletion re- 



spectively. It is interesting to note that f or c ase (a), 
where i^k = '^k = 1/a/2 and (/)k = 0, Eqs.(7[8) reduce 
to the ones for an ordinary condensate without SOC, 
where ng^/n ~ (na^)-^/^ and n|?/n ~ {na^Y^'^{T / for 
T < iii[17 . This fact allows a direct comparison among 
cases with or without SOC, as discussed later. 

Both the quantum and thermal depletion can be writ- 
ten as functions of a few dimensionless quantities. 



/n 



- Un)/U). (9) 
Un)/u.T/^), (10) 



where a = mil / (47r) is the averaged scattering length, 
and ^ = {2mnU)~^^^ is the healing length. Whereas the 
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values of Fq and Ft can be calculated numerically, it is 
useful to obtain an analytical expression of them using 
a qualitative analysis as follows. For quantum depletion 
in Eq.Q, the contribution from > /i to the integral 
is found to be rather small, similar to the case without 
S0C[17 . In the low energy regime, we further approx- 
imate (e+^q + e_^q)/2 by /i, its asymptotic value when 
q ^ 0. Then Eq.Q is reduced to 



Jo 



dep(e)/(2e), 



(11) 



where p{e) is the low energy DOS of quasi-particles de- 
termined by Eq.([5|. Similarly, for thermal depletion in 
Eq.([8|, the dominated contribution to the integral comes 
from < T, and Eq.([8| can be approximated by 



/iT / dep{e)/e^ 



(12) 



for low temperatures T < fi[22 . Eqs.(pT, 12) clearly 
demonstrate the important effect of p{e) on the conden- 
sate depletion. 

Using Eq.([6|, DOS of quasi-particles at low e for cases 
(a-c), which are represented by p'^^^(e), p^^\e) and p^^\e) 
respectively, can be obtained as shown in Table I. It 
is also useful to define an effective dimension according 
to the power-law behavior of p{e) . Considering a D- 
dimension system without SOC, the phonon-like quasi- 
particle (£^q = cq) gives rise to DOS as p{e) ~ e^~^. As a 
result, p^^\e)^ P^^\^) and p^^\e) correspond to effective 
dimensions D^ff = 3,5/2 and 2 respectively. The re- 
duced effective dimensions in cases (b) and (c) strongly 
indicates that the condensate depletion is significantly 
enhanced. 

Substituting the corresponding p{e) to Eqs.(ll, 12), 



we obtain the analytical expressions for condensate de- 
pletion in cases (a-c). As shown in the Tabel in Fig.l, 
both Fq and Ft take the form of scaling functions of a 
few microscopic variables. Two important results can be 
revealed from these forms. 

(i) Condensate depletion is largely enhanced in (b) and 
(c), as expected from the observation that D^ff is re- 
duced in these two cases. Compared with case (a), jn 
is enhanced by a factor of (/^^O^^^ (t>) and a fac- 
tor of for (c). Note that here K^i, = yjEf^j \i ^ 1. 
If Kx continuously increases, there is essentially no up- 
per limit for the enhancement. In current experiments, 
l/Ai^P"^, where A/, is the wavelength of the Ra- 
man field, and ~^/Al~5 — 10. Therefore, jn in- 
creases severalfold in case (b) and (c). For thermal deple- 
tion, the dependence of jn on temperature becomes 
~ (T//i)^/^ in case (b), in addition to the enhancement 
factor (j^x^^^^- Compared with the result ~ (T j p)^ in 
case (a), this means an even larger thermal depletion in 
case (b) at low temperatures T < p. 



(ii) The most significant effect of SOC occurs when 
considering the thermal depletion in case (c), which di- 
verges, i.e., ~ i^x^T J de/e. Therefore, the thermal 
depletion in three dimension becomes infinite in Bogoli- 
ubov theory at any finite temperature when both SOC and 
interaction become isotropic. This is a direct consequence 
of the reduction of D^ff to 2. 

It is worthwhile to point out that, while Bogoliubov 
theory is not expected to be accurate when the depletion 
becomes close to the total particle number, the above 
conclusion on the divergent thermal depletion for case (c) 
is consistent with the general argument that long-range 
order can not exist in a three dimensional system with 
two soft modes [23 . When the interaction can be written 
as U{n^ + n^)^, both the kinetic and interaction energy 
are constants on the sphere |k| = hcx- Two soft modes 
then appear along the tangent directions due to phase 
fluctuations in the long wavelength limit, whereas the 
density fluctuation is gapped and suppressed by repulsive 
interaction in this limit as ordinary condensates. This is 
similar to Rashba coupling in two dimension where one 
soft mode exists [24 . In our case, the two soft modes 
in phase fluctuation lead to divergent depletion and the 
absence of long-range order in case (c). 

To verify results (i) and (ii) , we perform numerical sim- 
ulations for the condensate depletion based on Eqs.([7|[8|. 
We first calculate the dependence of ng^/n and n|?/n on 
the anisotropy of SOC with a fixed value of tvx£,{= 8)- 
Fig. (2) shows how jn and /n evolve when the con- 
figuration of SOC changes from (a) to (b) by increasing 
Ky with Kz = ^ and then from (b) to (c) by increasing Kz 
with tvx = i^y One sees that jn increases from 0.15% 
in (a) to 0.5% in (b) and finally to 2.1% in (c). For ther- 
mal depletion, jn oX T / p = 0.5 increases from 0.09% 
in (a) to 1% in (b), and eventually becomes proportional 
io — — Kz/ i^x) when Kz approaches 

t^JX 9 as shown by 

the red dashed line in Fig. (2). We find that this loga- 
rithmical growth of jn exists for any finite value of 
T when SOC approaches the isotropic limit. Therefore, 
we conclude that jn is divergent at any finite tem- 
perature in (c). To further confirm the scaling forms 
in Fig.(l), we also calculate the dependence of conden- 
sate depletion in case (a-c) on dimensionless numbers in 
Eq.([lO|. As shown in the insets of Fig. (2), numerical 
results directly verify these scaling forms. 

We note that jn enters the logarithmic growth re- 
gion when 1 — Kzl ^ 10~^ for the parameters we used. 
While it might be challenging for current experiment to 
access this regime, the enhancement of thermal depletion 
is already visible far before entering this regime. As Fig. 2 
shows, n|?/n increases to 5% when Kzj^x = 0.98, which 
means n^/n has been enlarged by 50 times compared 
with case (a). Therefore, the enhanced depletion should 
be directly observable in experiments. 

Divergent depletion approaching case (c). An 
important question naturally arises that how the thermal 
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FIG. 2: (Color Online) ng^/n (upper panel) and /n at 
T I \i — 0.5 (lower panel) as functions of 1 — hiy/hix when — 
0(left column) and 1 — i^z/t^x when = /^y (right column), 
where U — U^i^ na^ — 10~^ and tix^ — 8. Fitting solid curves 
in orange are based on ^ 7^^^ for small Ky or a^^. Dashed 
red line is a linear fit ^ — ln(l — nz/i^x) when approaching 
isotropic limit. Inset of the upper panel is a In — In plot of 
In as a function of Kx^ for cases (a, b, c). Inset of the 
lower panel is a In — In plot of /n as a function of T//i for 
case (a,b). 



depletion evolves to infinity with increasing k^I^x from 
case (b) to (c), while keeping U = U^^^ fixed? Here we 
will demonstrate how an energy scale characterizing the 
anisotropy of SOC determines the depletion. 

For hiz ^ f^x = f^y and , the expansion of Eq. ^ along 
the z direction is well approximated by 



(o,o,g.) 




(13) 



where 7^ = 1 — /^^/a^^. Here 7^ reflects the suppression 
of the sound velocity along the z direction due to a fi- 
nite value of Hz- We define a characteristic momentum 
g* = 2^)1^ '^Kx by setting the contributions from and 



terms in Eq.(13) to be equal, i.e., ^zQ.'f' = 
g| corresponds to an energy scale = 2\^JzCK,x' 
£'* > /i > T in case (b), while ^ if case (c) is 
approached. Whereas both two terms in Eq.([l3| con- 
tribute to quasi-particle DOS, we use the following ap- 
proximation to obtain the analytic expressions for n|?/n 
in the crossover region from (b) to (c). For E < E^^ we 
ignore the contribution from q"^ term, and replace p(e) by 
p<(e) = jz ^^'^ p^^\e). For E > E*^ we ignore the contri- 



bution from term, and replace p{e) by p>(e) = p'^^^(e). 
Within this approximation, thermal depletion can be 
written as 

^ fiT r dep^{e)/e\ (14) 

JO 

n^^ ^ pT{f 'dep<(e)/e2+ / depy{e)/e^) (15) 
Jo Je: 



for E^ > T and E^ < T respectively. In another word, 
Eqs.( 14115 ) describe in the region with large and 
small value of 7^. Note that the above two equations 
become identical when E^ = T. To be more explicit, 
Eqs.( 14115 ) give rise to 



n; 



n) 



T 

I + /33 In 



T 



(16) 



,(17) 



where ft =1,2, 3 are numerical factors. 

Eq.(16) tells one that, starting from case (b), ther- 
mal depletion increases as a function of 7^ for small 
values of Kz- This behavior of thermal depletion is di- 
rectly verified by the numerical results in Fig. (2), as 
demonstrated by the fitting curves in orange. For a large 
value of Kz ^ ^ic, Eq.([T7|) clearly shows that, at any 



given temperature, becomes logarithmically diver- 
gent ~ — ln£^* ~ —hi{l — Kz/ i^x) when Kz approaches 
Therefore, Eqs.([l6j [l7| can be used to qualitatively 
describe the thermal depletion in the crossover region 
from (b) to (c). Similar discussions apply to the conden- 
sate depletion in other crossover regions [21]. 

In conclusion, we have pointed out that SOC signif- 
icantly enhances condensate depletion of bosons, and 
completely destruct the condensate when both SOC and 
interaction approach the isotropic limit. We also show 
that the dependance of condensate depletion on micro- 
scopic variables is distinct for different types of SOC. 
These results are expected to be useful for future ex- 
perimental exploration of spin-orbit coupled condensates, 
and will stimulate more studies on strongly depleted con- 
densates for which correlation effects are important. 
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SUPPLEMENTARY MATERIAL 

In this supplementary material, we present the re- 
sults for mean- field ground states, and analytical expres- 
sions of both the quantum and thermal depletion for the 
crossover regions from (a) to (b), and from (b) to (c). 



Mean-field ground state 

The kinetic energy shown in Eq.(2) of the main text 
leads to different manifolds of single-particle ground 
state, as listed below. 

{1) > i^y^ i^z ' two points ±ko, with kg = /^x^x 

(2) K/j; K/y ^ K/2; ' circle ~\~ ky ^x") ~~ 

(3) f<ix = i^y = i^z ' ^ sphere |k| = k^x 

For (1) and (2,3), the single-particle ground states are 
doubly and infinitely degenerate respectively. The pres- 
ence of interaction may lift these degeneracy. Consider a 
plane- wave condensate |P)k = <^k^|0) and a stripe con- 
densate l^*)!;: = ( ^^^~^ )^|0), their interaction energies 
can be written as 

Ep{k) = {P\Hi^t\P) = ^{U + 2{Un - U)ulv^^ 

i?5(k) = {S\Hi^t\S) = ^({U + Uu) - 2(C/n - U)ulvi 

The expressions for and have been given in the 
main text and '^k + = 
determined as follows. 

(1) As «|ko = «Lo 
^^^(ko), we see that the ground state is |P)±ko fo^" U > 
U^l, or |6')ko for U < U^i. 

(2) The spontaneous symmetry breaking can occur 
along any direction in the x-y plane. The ground state is 
|P)k^ and |5')k^ for U > U^^ and U < U^^ respectively, 
where k^ is an arbitrary rotation of ko about the z axis. 

(3) Despite the high degeneracy of kinetic energy for 
this case, expressions of £^p(k) and Es(k) tell one that 
interaction energy is minimized when the value of u^v^. 
is either 1 or 0, i.e., when k = k^ or k = ±ki with 
ki = f^z^z- Compare £^p(k0), ^^^(k^), £^p(±ki) and 
^S'(ki), we find the ground state is |P)k0 for U > U^i^ 
and |P)±ki for U < U^^. It is worthwhile to mention 
that though |5')ki is degenerate with |P)k0 for U > U^^, 
it is physically unachievable, as spin- up and spin-down 
particles cannot be mixed with each other k^ = ky = 0. 

The results for (1-3) are summarized in the following 
table. 



1. The ground state for (1-3) is 
= 1/2, compare £^p(±ko) and 
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(1) 


(2) 


(3) 


U>Un 








U <Uu 









fiTj^Uep^{e)/e' and ^ pT{jfy dep^{e)/e' ^ 

jl^ dep^{e) / e^) for > T and < T respectively, 
which give rise to 

nf?/n = a;7;M^^')"^(W, (23) 
nf?/n = (l - ^^3(^,7^)) ,(24) 

where a-^^^ 2 3 numerical factors. 



TABLE I: Mean-Field Ground state in 3D. 

Condensate depletion for the crossover 

(A) From (a) to (h) 



For = {) and an arbitrary ny/^x ^ [O,!], Eq.(5)(in 
the main text) at low momentum along the y direction 
can be expanded as 



E^ 



(0,gy,0) 



where 



_ 1 



(19) 



H (1 - x^) 



2\2 



n ^2 



x\l-x^). (20) 



For small tiy/tix <^ 1^ Eq 18 is dominated by the lin- 



ear mode, and the sound velocity along the y direction is 

]^ /2 ' 

reduced by a factor of . li Ky/ ^ 1^ Eq : 



18 



is dom- 
inated by the quadratic mode. As p <^ E^, the crossover 
from the linear mode dominated to quadratic mode dom- 
inated regimes occurs when I Kx 1, and 
f{x) 1. In this limit, by setting jyQy'^ = g'*^/(4/^^), we 
define a momentum scale to characterize the crossover, 



El 



o 1/2 



E^ 



(o,<?;,o) 



and the corresponding energy scale is 
= 2y/2^ycnx. 



Similar to the discussions in the main text, we ignore 
the term in Eq.(18) for £^ < in calculating the 
analytic form of quasi-particle DOS, and obtain p<(e) = 
^y^^'^ p^^\e)' In contrast, for E > Ey, we ignore Qy and 
obtain p>(e) = p^^\e). Within this approximation, quan- 
tum depletion can be written as ng^ ~ p dep^{e) / e 

and ng- « /x(/o^" dep<{e)/e + dep>(e)/e) for /x < E; 
and p > Ey respectively. To be explicit, we obtain 

ng^/n = ai7| (na^) ^ , (21) 
n^^yn = a2{^xOHna^y^ (l - as{E;/ p)l^ , (22) 

for these two cases, where a^^i^2,3 are numerical factors. 
When Ey decreases down to zero at Kjy Kjx 1 

Eq.rt22k 

becomes the expression for Rashba coupling, i.e, jn ~ 
(^a^O ^ (na^) 2 , which is consistent with a numerical result 
found in Ref.[T5j. 

The same discussions apply to thermal depletion 
for this crossover, which can be written as ~ 



(B) From (b) to (c) 

For Hvy Hvx 

and an arbitrary tZz/f^x ^ [O7 1]^ Eq.(5)(in 
the main text) at low momentum along z-direction can 
be expanded as 



E^ 



m 



m..) = J^(<lh^ + ^9i-)] (25) 



4«r 



where 



_ /t| n{U-UH) . 



P 



(26) 



(27) 



Similar to previous discussions for the crossover from 

(a) to (b), here by setting 7^^^^ = g*^/(4Ai:^), we define 

a characteristic momentum for the crossover from (b) to 
1/2 

(c), = 27^' Kx, and the corresponding energy scale 
El = 2^/2^,CKx. 

Following the same procedure in previous discussions, 
quantum depletion can be written as 

nr/n = /^;7z"*(KxC)*(na3)i (28) 
n^/n = /3^(/t,e)(na3)5 (l - /3^(E*/m)*) (29) 

for E^ > p and E^ < p respectively, where Pl^-^ 2 3 
are numerical factors. With decreasing the scaling 



form of the quantum depletion involves from Eq.(28) to 
Eq.([29|) dit El = p and eventually becomes the one for 
case (c) with isotropic SOC and interaction. 

At finite temperatures, the thermal depletion can be 
written as Eq.(16) and (17) in the main text for E^ > T 
and El <T respectively. At El = T, the two expressions 
are identical. With further increasing El down to zero, 
thermal depletion eventually becomes divergent, in the 
form of — ln(l — hCz/f^x) if /7 = U^i as shown in Fig. 2 in 
the main text. 



